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Oxygen vacancies ( V Q ) effects on magnetic ordering in Eu 0 .5Ba 0 .5TiO 3 _5 (EBT0 3 _s) thin films have been 
investigated using a combination of experimental measurements and first-principles density-functional 
calculations. Two kinds of EBT0 3 _§ thin films with different oxygen deficiency have been fabricated. A 
nuclear resonance backscattering spectrometry technique has been used to quantitatively measure contents 
of the V 0 . Euo.5Ba 0 .5Ti0 3 ceramics have been known to exhibit ferroelectric (FE) and G-type 
antiferromagnetic (AFM) properties. While, a ferromagnetic (FM) behavior with a Curie temperature of 
1.85 K has been found in the EBT0 3 _§ thin films. Spin-polarized Ti 3+ ions, which originated from the Vq, 
has been proven to mediate a FM coupling between the local Eu 4f spins and were believed to be responsible 
for the great change of the magnetic ordering. Considering the easy formation of Vq, our work opens up a 
new avenue for achieving co-existence of FM and FE orders in oxide materials. 



Multiferroics, or materials that simultaneously possess two or more ferroic orders [ferroelectric, (anti-) 
ferromagnetic, and ferroelastic] , have recently returned to the forefront of materials research due to their 
rich physical properties and potential applications in data storage, sensors, and spintronics 1 ' 2 . From an 
application point of view, multiferroics with ferromagnetic-ferroelectric (FM-FE) orders are more attractive and 
highly desired 3 . However, very few exist in nature due to the intrinsic contradiction in existence between the FM 
and FE ordering within a single phase 4 . On the other hand, antiferromagnetic-ferroelectric (AFM-FE) materials 
are more commonly found, such as rare-earth manganites (RMn0 3 ) and the well known BiFe0 3 511 . It is 
fundamentally interesting and technologically important to develop FM-FE materials by changing magnetic 
ordering of the AFM-FE materials. An additional driving force and a deep understanding of physical phenom- 
enon underlying the transformation are needed to realize this purpose. 

In oxides, oxygen vacancies (V 0 ) have been approved to be intrinsic defects and are believed to have a critical 
impact on their properties 12 . Coey et al reported that V Q play a key role in obtaining room temperature FM in 
Hf0 2 , ZnO, Ti0 2 , and other non-magnetic oxide systems 13 " 16 . It has also been proven that V Q are useful in 
increasing Curie temperature and enhancing magnetic moment of EuO thin films 17 " 20 . Furthermore, V Q are 
known to induce a room temperature ferroelectricity in SrTi0 3 thin films 21 " 23 . Despite the fact that the V Q are 
effective to manipulating the magnetic and ferroelectric properties of oxides, little attention has been focused 
upon the V Q effects on multiferroicity in single-phase materials. A natural question to ask is: are the V Q able to 
change the magnetic ordering from AFM to FM in multiferroic materials? 

The beginning of the present work is to choose a suitable material to verify the above idea. In recent years, 
divalent europium oxidation materials showed attractive functionalities, rendering them subjects of intensive 
studies. Compared to the 0.05 fi B /¥e in BiFe0 3 24 , rare-earth Eu 2+ ion has seven unpaired and localized 4f 
electrons, which resulted in a large magnetization of 7 /n B /Eu and sometimes is coupled with electrical properties. 
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For example, Rushchanskii et at confirmed that Eu 0 .5Ba 0 .5TiO 3 
(EBTO) ceramics exhibit FE (Curie temperature T c ~ 213 K) and 
G-type AFM (Neer temperature T N ~ 1.9 K) properties 25 . And coup- 
ling between magnetism and dielectric properties was predicted in 
Eui_ x Ba x Ti0 3 materials 26 . Furthermore, as we discussed earlier, the 
V Q have been proven to be useful in manipulating the magnetic 
orders of divalent Eu ions in EuO thin films. Therefore, EBTO is a 
potential material for the present study, based on the nature of a large 
magnetization of 7 /n B /Eu and the AFM-FE orders. In addition, if we 
can alter magnetism of EBTO thin films by doping V Q and at the 
same time still preserve good ferroelectricity, it will certainly make 
EBTO itself appealing in fundamental research and practical appli- 
cations. In the present work, we experimentally confirmed that the 
magnetic ordering in oxygen- deficient Euo.5Bao. 5 Ti0 3 _§ (EBT0 3 _§) 
thin films has been transformed from AFM to FM. First-principle 
calculations indicated that spin-polarized Ti 3+ ions, which originated 
from the V 0 , mediated a FM coupling between the local Eu 4f spins. 
Considering the easy formation of V 0 , our work presented an effec- 
tive technique to achieve co- existence of FM and FE orders in oxide 
materials. 

Results 

The EBT0 3 _ 5 thin films were fabricated on (001) oriented SrTi0 3 
(STO) and Nb-doped SrTi0 3 (Nb-STO) substrates by pulsed laser 
deposition (PLD). Following the deposition, parts of the thin films 
were post- annealed under a reducing atmosphere to increase amount 
of the V Q . Typical X-ray diffraction (XRD) 0-20 scans confirmed 
that the EBT0 3 _§ is preferentially oriented along the c axis (not 
shown). In order to reveal the microstructure of the EBTO 3 _§ thin 
films, cross-sectional transmission electron microcopy (TEM) mea- 
surements were conducted. A low magnification bright- field TEM 
image (Fig. la) shows that the interface is clean without any obvious 
interface reaction or intermixing. In the view area, the film is con- 
tinuous without any obvious grain boundaries, confirming the 2D 
growth of the film. The corresponding selected area diffraction 
(SAD) image taken from the interface area is shown as inset of 
Fig. la. The distinguished diffraction dots from the film and substrate 
indicate the high quality epitaxial growth. The out-of-plane and 
in-plane orientation relationships have been determined to be 
(001) E bto|| (OOI)sto and [110] EBT o||[110]sto> respectively. A high 
resolution TEM (HRTEM) image (Fig. lb) confirmed the excellent 
heteroepitaxial relation with an atomically sharp interface. Similarly 
no obvious misfit dislocations have been found along the interface, 
which is consistent with the SAD observation. 

To quantitatively analyze stoichiometry and oxygen concentration 
of the EBT0 3 _s thin films, a nuclear resonance backscattering spec- 
trometry (NRBS) technique was used. The advantage of NRBS over 
traditional Rutherford backscattering spectrometry (RBS) is its sens- 
itivity to selectively measure oxygen in the EBT0 3 _s film on an 
oxygen-rich substrate, such as STO in our experiment 27,28 . The 
NRBS was performed using a 3.043 MeV 4 He + analyzing beam 
and the backscattering particles were detected by a surface barrier 
silicon detector located at 167° from the beam direction. At such a 
beam energy, the He scattering from the heavier elements (Eu, Ba, 
and Ti) in the film is still Rutherford, therefore the Eu:Ba:Ti ratios 




Figure 1 | TEM images of the EBT0 3 _s thin film on a STO substrate: 
(a) low magnification cross- sectional TEM; (b) high resolution TEM 
(HRTEM). The inset shows the corresponding selected-area diffraction 
(SAD) pattern. 

can be determined reasonably well by fitting the experimental spectra 
data with the commercial Rutherford Universal Manipulation 
Program (RUMP) software 29 . Within the uncertainty of the measure- 
ments (—5%), we determined that our EBT0 3 _§ thin films do have 
an Eu:Ba:Ti ratio of 1 : 1 : 2. 

After the atomic ratios of cations were determined, a nuclear 
scattering and reaction simulation package (SIMNRA) was used to 
fit the nuclear resonant oxygen scattering spectra (not shown) 30 . To 
minimize uncertainties related to the nuclear scattering cross section 
and the incident beam energy, a bare STO substrate as a standard 
reference was measured and analyzed all together along with our 



Table 1 | The measured concentration of europium and oxygen from the NRBS. The content of Vo(x), thevalueof 5 (5 = 3x), and the ideal 
and real content of Ti 3+ are also shown. The ideal content of Ti 3+ is estimated from 5 with a relationship of Ti 3+ /Ti:5 = 2:1. The real values 
are calculated by fitting peaks of XPS measurements 

Oxygen vacancies 

Eu(at./cm 2 ) 0(at./cm 2 ) x 5 Ideal Ti 3+ /Ti Real Ti 3+ /Ti 

As-deposited thin film 8.7xl0 16 5.01 Xl0 17 4.0% 0.120 24.0% 27.1% 

Annealed thin film 9.8xl0 16 5.58 X10 17 5.1% 0.153 30.6% 33.4% 
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Figure 2 | The open circles represent peak fittings on these graphs, which 
show the Ti 2p core-level photoemission spectra of (a) as-deposited and 
(b) annealed EBT0 3 _s thin films. 
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Figure 3 | Temperature dependence of magnetization curves under ZFC 
and FC conditions for the (a) as-deposited and (b) annealed EBT0 3 _g 
films. The insets show the derivative of the magnetization (the solid 
circles) and the reciprocal susceptibility (the open circles) with respect to 
the temperature (obtained from the FC curves). The blue lines on the insets 
are the Curie-Weiss law fittings of the reciprocal susceptibilities. 



EBT0 3 _s thin films. Table 1 shows the measured concentration of 
europium and oxygen in the as-deposited and annealed EBT0 3 _5 
thin films. From the Eu concentration, the ideal oxygen concentra- 
tion in the stoichiometric Eu 0 .5Ba 0 .5TiO 3 thin films is estimated to be 
5.22 X 10 17 and 5.88 X 10 17 at./cm 2 for as-deposited and annealed 
thin films respectively. By comparing the measured and the ideal 
oxygen concentration, we have calculated that there are 4.0% oxygen 
deficiency (5 ~ 0.120) in the as-deposited film and 5.1% oxygen 
deficiency (5 ~ 0.153) in the annealed film. In other words, the V Q 
were introduced into the EBT0 3 _§ thin films and the content of the 
V Q were increased by post- annealing under a reducing atmosphere. 

Because of the charge compensation, the existence of V Q should 
induce changes in the valance states of cations. To quantify these 
changes, we investigated the valance states of the EBT0 3 _5 thin films 
by using X-ray photoemission spectroscopy (XPS). We found that 
the valence states of Eu and Ba ions are remained in divalent (not 
shown), while those of Ti ions have been changed. Because V Q act as 
rc-type dopants, the Ti 3+ and Ti 4+ formal valences should coexist in 
the EBT0 3 _§ thin films. Figure 2 shows Ti 2p core-level photoemis- 
sion spectrum of (a) the as-deposited and (b) the annealed EBT0 3 _ 5 
thin films with peak fittings. Compared with previous results, the 
peaks at a binding energy around 459.42 eV and 465.45 eV can be 
assigned as Ti 4+ , with peaks around 457.81 eV and 463.53 eV cor- 
responding to Ti 3+ 31 ' 32 . By comparing the peak area, the content of 
Ti 3+ was calculated to be 27.1% and 33.4% in the as-deposited and 
annealed thin films respectively. On the other hand, by assuming two 
free carries for each Vq 33 , the ratio between the content of Ti 3+ and 
that of V 0 (5) will be 2 : 1. Therefore, the ideal content of Ti 3+ for the 
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Figure 4 | Magnetic field dependent magnetization curves for (a) as- 
deposited and (b) annealed EBT0 3 _§ thin films at various temperature. 

The insets show magnetic susceptibility curves. 
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Figure 5 | Spin-density and local density of states (DOS) in the EBT0 3 _ 1/8 : (a) and (c) for FM ground state; (b) and (d) for AFM state. Here, yellow 
color is for up-spin and blue color for down-spin. The red vertical dash line is the Fermi level. Til, Ti2, and Ti3 are the nearest (just below the V 0 ), 
next-nearest (at the same ab-pfane of Til), and next-next-nearest (at the same afr-plane of Til) Ti ions with respect to the V Q . The Eu ions are almost 
similar to each other. The development of spin-polarized d z 2 orbital at Ti 3+ is obvious in the FM ground state. 



as-deposited and annealed thin films should be 24.0% and 30.6% 
respectively (see table 1). These results are well fitted with the calcu- 
lated values and further confirmed that V Q were doped into the thin 
films. In addition, the fact that the change of the valence states only 
occurred in Ti cations has been verified by first-principles density- 
functional calculations and will be discussed in more details later. 

To investigate V Q effects on magnetic properties of the EBT0 3 _§ 
thin films, temperature and magnetic-field dependences of magnet- 
ization were investigated using a superconducting quantum interface 
device magnetometer (SQUID). Figure 3 shows the temperature 
dependent magnetization curves. The measurements were per- 
formed under zero-field-cooled (ZFC) and field-cooled (FC) condi- 
tions with an external magnetic field of 100 Oe applied parallel to the 
films surface. The EBTO ceramics antiferromagnetically ordered at 
1.9 K 25 . However, the magnetization on both ZFC and FC curves 
increased monotonically with decreasing temperature until 0.9 K 
and tend to be saturated at the lowest temperature, which is a typical 
FM behavior. The derivative of the magnetization and the reciprocal 
susceptibility as a function of temperature are shown as insets of 
Fig. 3a and 3b. It has been found that the derivative of the magnet- 
ization shows a sharp valley, in other words the magnetizations 
increase quickly, at around 1.85 K. In addition, the Curie- Weiss 
law fitting of the reciprocal susceptibilities at high temperatures 



intersects the temperature axis at 1.85 K. Both of these results then 
confirmed the FM ordering in the EBTO3-5 thin films with a Curie 
temperature of 1.85 K 34 . Figure 4 shows the field dependent magnet- 
ization curves at 1, 1.5, and 5 K. It has been found that, at 1 and 1.5 K, 
the magnetization increased sharply at low fields and then 
approached saturation quickly. The large susceptibility at low fields 
(shown as insets of Fig. 4a and 4b) is also an evidence for the FM 
ordering below 1.85 K. The saturation magnetizations, measured for 
the as-deposited and annealed thin films at 1 K, are around 6.75 
Eu and 6.85 fi B /Eu respectively. These values are close to the ideal 
magnetic moment of the Eu 2+ ions (7 ^b/Eu). It should be noticed 
that the susceptibility at low fields and the value of the magnetization 
in the annealed thin films are larger than those in the as-deposited 
thin films, which obviously originates from the enhanced content of 
V Q . In order to further prove the FM ordering of the EBT0 3 _s thin 
films below 1.85 K, magnetic hysteresis loops were measured at 
0.5 K (Supplementary Figure SI). Hysteretic behavior is observed 
with a coercivity of 30 Oe. The combination of these results proved 
that the EBT0 3 _§ thin films become a ferromagnet with a Curie 
temperature of 1.85 K, which is different from bulk EBTO - G-type 
antiferromagnet 25 . 

Many theoretical studies of vacancy-induced magnetism in 
non-magnetic SrTi0 3 have been performed using ab initio 
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Table 2 | Energy differences between AFM and FM states in EBT03_ 1 74 with various atomic orderings of Eu and Ba ions. In all the cases, FM 
state is energetically favorable. In addition, A-type atomic arrangement of Eu and Ba ions is the ground state 

A-type atomic ordering C-type atomic ordering G-type atomic ordering 



AE = E(AFM) - E(FM) (meV/Eu) 0.6 0.2 0.5 



calculations 35 ' 36 . To investigate the fundamental physics underlying 
the origin of the FM ordering in the EBT0 3 _s thin films, we calcu- 
lated total energies of EBT0 3 _ 1/8 with both FM and AFM orderings 
from first-principles density- functional theory. FM ordering is found 
to be 13.1 meV/Eu favorable. Figure 5 shows the spin-density and 
local density of states (DOS) in the EBT0 3 _ 1/8 : (a) and (c) for FM 
ground state; (b) and (d) for AFM state. As shown in DOS of local Ti 
ions, the occupancy of 3d orbital indicates the appearance of Ti 3+ , in 
good agreement with the XPS measurements. Furthermore, the exist- 
ence of V 0 does not change the valence state of Eu, since almost fully 
occupied 4f 7 states are found, which is similar to the pristine EuTi0 3 
and confirmed by the SQUID and XPS measurements 37 . On the other 
hand, as shown in Fig. 5 (a) and (c), the development of large local 
spin moment at Ti sites makes the FM ground state differ signifi- 
cantly from the AFM state. The local magnetic moment at Ti 3+ ions 
right above or under V Q is around 0.24 fi B and the d z 2 dominates. 
While for other Ti 3+ ions, the local orbital is either d xy or d xz /d yz . 
Magnetic moment at Eu site does not change from the pristine case 
with a local magnetic moment of around 7 fi B . Judging from Fig. 5, it 
is the appearance of spin-polarized Ti 3+ that mediates the FM coup- 
ling between the localized Eu 4f spins. It should be noticed that with 
increase of the content of V Q , e.g. Eu 0 .5Ba 0 .5TiO 3 _ 1 /4, the FM order- 
ing is still found to be energetically favorable than the AFM orderings 
(see Table 2). These results verified that the EBT0 3 _§ thin films do 
show FM ordering, which originated from the ordering of the spin of 
Eu 4f electrons mediated by the spin-polarized Ti 3+ ions. 

This article is mainly concerned about the magnetic ordering, but 
we also investigated the FE properties of the annealed EBT0 3 _s thin 
films by measuring dielectric constant as a function of temperature 
and room-temperature hysteresis loop using piezoresponse force 
microscopy (PFM). It's surprised to find out that a peak in dielectric 
constant versus temperature curves (shown as Fig. 6) is clearly seen 
above room temperature. Moreover, a room-temperature hysteresis 
loop measured by PFM was shown as inset of Fig. 6. These results 
suggest that the FE Curie temperature of the annealed EBT0 3 _ 5 thin 
films is above room temperature, which is significantly higher than 
that of bulk EBTO (—213 K). The investigation on mechanism of the 
enhancement on the FE Curie temperature is presently in progress. 
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Figure 6 | Temperature dependence of dielectric constant for the 
annealed EBT0 3 _s thin film. The inset shows room-temperature 
hysteresis loop measured by PFM. 



Discussion 

In summary, two kinds of EBT0 3 _§ thin films with different content 
of V Q were fabricated. By using NRBS, SQUID, and first-principles 
calculations, we demonstrated that the magnetic ordering of oxygen- 
deficient EBT0 3 _5 thin films has been changed from AFM to FM. 
The transformation originated from the V Q induced Ti 3d electrons, 
which mediated the FM coupling between the local Eu 4f spins. The 
results are significant as they proved that V Q are effective to manip- 
ulating magnetic ordering in multiferroic materials. Considering the 
easy formation of V 0 , the present work presents a methodology to 
enhance multiferroicity in the EBT0 3 _s thin films and this method 
holds great promise for other oxide materials. 

Methods 

The EBT0 3 _ 5 thin films were fabricated on (001) oriented SrTi0 3 (STO) and Nb- 
doped SrTi0 3 (Nb-STO) substrates by pulsed laser deposition (PLD) using a pulsed 
excimer laser (Lambda Physik, 248 nm, 3 Hz, 2 J/cm 2 ). A high-density EBTO 
ceramic pellet was used as the target. The details of the preparation of EBTO ceramics 
can be found in elsewhere 25 . Deposition temperature was 700°C and oxygen pressure 
was 1 X 10" 4 Pa with the purpose of doping V Q . Following the deposition, parts of the 
thin films were annealed at 1000°C under a flowing gas of 95 vol% Ar + 5 vol% H 2 for 
10 hours to increase amount of the V Q . The film thickness, revealed by cross- sectional 
transmission electron microcopy (TEM), was 100-150 nm. 

The crystal structures were characterized by X-ray diffraction (XRD, Rigaku K/ 
Max) and TEM (FEI Tecnai F20 analytical microscope). Nuclear resonance back- 
scattering spectrometry (NRBS) was performed on Los Alamos National Laboratory. 
A He + beam energy of 3.043 MeV was used to quantitatively analyze stoichiometry 
and oxygen concentration in the EBT0 3 _§ thin films. The valence states were 
investigated by X-ray photoemission spectroscopy (XPS) at PHI5000 VersaProbe. 

Magnetic measurements were performed on EBT0 3 _g/STO thin films using a 
superconducting quantum interface device magnetometer (SQUID) equipped with a 
He 3 insert (Quantum Design, MPMS-XL). The electrical properties were measured 
using a Pt/EBT0 3 _§/Nb-STO heterostructure. The dielectric constants were inves- 
tigated using an Agilent 4294A Impedence Analyzer. The measurements were per- 
formed at selected temperatures in a Linkam Scientific Instruments HFS600E-PB4 
system. Room-temperature piezoresponse force microscopy (PFM) was measured 
using PFM mode of Asylum Research MFP-3D-SA atomic force microcopy. 

Our ab initio calculations are performed using the accurate full-potential pro- 
jector-augmented wave (PAW) method, as implemented in the Vienna ab initio 
simulation package (VASP) 3839 . They are based on density- functional theory with the 
generalized gradient approximation (GGA) in the form proposed by Perdew, Burke, 
and Ernzerhof (PBE) 40 . The on-site Coulomb interaction is included in the GGA + U 
approach with effective U = 4 eV for Eu 4/orbitals 41 . A plane-wave cutoff of 600 eV is 
used throughout and the convergence criteria for energy is 10" 6 eV. PAW potentials 
are used to describe the electron-ion interaction with 17 valence electrons for Eu 
(4f 5s 2 5p 6 6s 2 ), 10 for Ba (5s 2 5p 6 6s 2 ), 10 for Ti (3p 6 3d 2 4s 2 ), and 6 for O (2s 2 2p 4 ). In our 
calculations, ions are relaxed toward equilibrium positions until the Hellman- 
Feynman forces are less than 1 meV/A. In addition, lattice constants are optimized 
until the stress is less than 0.001 Pa. The fully optimized lattice constant of 
Eu 0 .5Ba 0 .5TiO 3 stems from an initial cubic lattice, from which the atoms shift a little 
from their high symmetric positions. Total energies are calculated and compared 
between different magnetic and A-site atomic arrangements. 
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